Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
黒澤 亮平; 岡田 翔太; 坂井 章浩; 仲田 久和; 天澤 弘也
JAEA-Data/Code 2015-005, 82 Pages, 2015/06
研究施設等廃棄物のうち国内の研究用原子炉施設の解体に伴って発生する廃棄物の放射能インベントリ評価に使用する原子炉内及び周辺の中性子束の計算を汎用的かつ簡易に行える中性子束計算ツールを整備した。本ツールは、巨視的実効断面積の計算、中性子束の計算及び放射化イベントリの評価するために算出した中性子束の編集を行う計算プログラムから構成される。本報告書は、研究用原子炉施設における中性子束の評価方法の概要、中性子束計算ツールの全体構成と機能、各計算プログラムの入出力データと試計算をまとめたものである。
北端 琢也; 清田 史功; 白鳥 芳武; 井口 幸弘; 松井 祐二; 佐藤 裕之
JNC TN3410 2000-014, 43 Pages, 2000/09
新型転換炉ふげん発電所は核燃料サイクル開発機構法にもとづき、平成15年度までに運転を停止することになっており、現在、廃止にともなう措置に関する技術の開発及びこれに必要な研究(以下「廃止措置技術開発」という)を実施している。この廃止措置技術開発を計画・実施するにあたり、「ふげん」を国内外に開かれた技術開発の場として十分に活用するとともに、当該技術開発で得られる成果を有効に活用することを目的として、サイクル機構外の有識者で構成される「ふげん廃止措置技術専門委員会」を平成11年12月に設置し、平成11年12月14日に第1回委員会が開催された。平成12年度も引き続き設置され、平成12年8月28日に第2回委員会が開催された。本書は、第2回ふげん廃止措置技術専門委員会において配布された資料集であり、「ふげん」廃止措置への取り組み状況、生体遮へい体コンクリートの放射化量評価、コンクリート中のトリチウム濃度測定方法の検討、「ふげん」の廃止措置技術開発の進め方についてまとめたものである。また、併せて、当該委員会において参考として報告された系統化学除染の実施状況についても記載した。
庄野 彰; 佐藤 若英*; 栗原 国寿
JNC TN9400 2000-037, 87 Pages, 2000/03
水冷却型増殖炉の核特性の特徴を把握するために、公開文献に基づいて沸騰軽水を冷却材とした炉心、それをベースに冷却材を非沸騰重水及び非沸騰軽水に置換した炉心、ならびに高速炉の代表炉心として大型Na冷却MOX燃料炉心の仕様を設定し、基本的な核特性の比較検討を実施した。高速炉用非均質セル計算コードSLAROMと軽水炉解析に汎用的に用いられているSRACの解析結果を比較し、転換比・中性子スペクトル・エネルギー領域別反応割合・1群断面積等について、両コード間の差は小さく、高速炉核特性評価システムが水冷却型増殖炉の基本的核特性の検討に適用可能であることを確認した。SLAROMコードを用いて上記4種類の炉心の核的パラメータ及び増殖性に見られる相違を考察した。冷却材の変更によって、中性子スペクトル・値・主要エネルギー領域等が変化する傾向を把握した。水冷却型増殖炉では、冷却材中に存在する水素(または重水素)の影響で低エネルギー成分の中性子束が高速炉に比べて増大し、その結果、主要エネルギー領域が低エネルギー側に移行し、核分裂性核種の値が低下するが、MOX燃料を稠密に配置して増殖性を担保する設計であることが理解できた。Pu同位体組成が増殖特性に及ぼす影響をSRACコードの燃焼解析機能を用いて検討し、Pu-240含有率が大きい組成のPuを装荷した場合には転換比が大きく算定される傾向にあることを定量的に評価した。臨界性及び反応率比に対する感度解析により、沸騰軽水冷却増殖炉では、高速炉に比べて、1KeV以下のエネルギー領域における感度が増加することがわかった。断面積不確かさに起因する核設計予測精度評価については、現在核データセンターで共分散データの見直し中であるため、現状の共分散データを用いて暫定解析結果を得た。見直し後の最新共分散データを用いた内部転換比予測精度評価が今後の課題である。
竹村 守雄*
JNC TJ9450 2000-002, 112 Pages, 2000/03
本報告書は、1986年に開始されたJASPER(Japanese-American Shielding Program for Experimental Research)計画の中で企画された計8個の一連の実験の中の6番目の実験として、1992年3月初めから約2ケ月間かけて米国オークリッジ国立研究所(ORNL)にて実施されたギャップストリーミング実験の測定値、実験体系の寸法・組成データ、実験状況、測定システム等の情報を一冊にまとめたものである。作成にあたっては、ORNLから発行されたギャップストリーミング実験結果報告書の内容を基本とし、それに現地派遣員より報告された情報を補足した。ギャップストリーミング実験は、高速炉原子炉容器上部の中性子束がインクロージャシステムに存在する空隙部を通じての中性子ストリーミングに起因していることから、空隙部の幅、オフセット間隔と中性子ストリーミングの大きさの相関を実験的に把握すると共に、このストリーミング評価に適用される解析手法の精度の検証と向上のためのデータベースを提供するため企画された。ORNL遮蔽実験施設TSF(Tower Shielding Facility)の実験炉TSR-II(Tower Shielding Reactor-II)のコリメータ後方に鉄ライナー付きコンクリート供試体が設置された。コリメータと供試体との間にスペクトルモディファイアを設ける体系と設けない体系とが構成され、軟らかいエネルギースペクトル(実機を模擬)と硬いスペクトルの中性子が供試体に入射された。コンクリート供試体には中央に貫通部があり、コンクリート製シリンダー、円筒状スリーブを挿入組合わせることにより、円環状ギャップの幅とオフセット位置を変えられる構造となっており、供試体の後方において各種検出器を用いた中性子測定が行われた。広いエネルギー範囲にわたる中性子が、ビーム軸に直交する方向での分布及びビーム軸上において、ほぼ全体系において測定された。硬い入射スペクトルのものについてはさらに、高速中性子のビーム軸上でのエネルギースペクトルの測定及びビーム軸に直交する方向の分布の測定が行われた。なお、JASPER計画は日米両国の液体金属冷却炉(LMR)の現行設計内容の進展に資することを目的として、米国エネルギー省(USDOE)と動力炉・核燃料開発事業団(現核燃
中島 宏; 益村 朋美*; 田中 進; 坂本 幸夫; 高田 弘; 田中 俊一; 中根 佳弘; 明午 伸一郎; 中村 尚司*; 黒沢 忠弘*; et al.
Journal of Nuclear Science and Technology, 37(Suppl.1), p.192 - 196, 2000/03
日本原子力研究所TIARAの第2軽イオン室において、67MeV陽子のp-Cu反応から生じる中性子のターゲット周囲における角度・エネルギー分布を放射化法により測定した。また、中性子ターゲット室内における速・熱中性子束の空間分布を核分裂計数管、TLD及び金の放射化法によりそれぞれ測定した。今回は、この実験について、MCNP-4B等のモンテカルロ計算コードによる解析を行い、その計算精度を検証した。その結果、中性子線源であるターゲット構造体周囲の速中性子束については、後方角を除いて、計算値は角度・エネルギー分布の測定値と50%以内で一致した。室内の速中性子束分布については、線源から離れ、散乱線の寄与が大きくなるにつれて、計算値が過小評価する傾向にあるが、50%以内で測定値と一致した。一方、熱中性子束分布については、計算値は全体的にファクター2程度測定値を過大評価した。
中島 宏; 益村 朋美*; 田中 進; 坂本 幸夫; 田中 俊一; 中根 佳弘; 明午 伸一郎; 中村 尚司*; 黒沢 忠弘*; 平山 英夫*; et al.
Journal of Nuclear Science and Technology, 37(Suppl.1), p.197 - 201, 2000/03
日本原子力研究所TIARAの第2軽イオン室に続く迷路状通路において、67MeV陽子のp-Cu反応による中性子を用いて、中性子エネルギースペクトルと熱中性子束及び線量当量空間分布を測定した。この実験について、モンテカルロ計算手法及び簡易計算手法により解析し、その精度を検証した。その結果、線量当量分布に関して、第1脚でモンテカルロ法による計算値がファクター2程度過小評価し、第3脚では過大評価となったが、第2脚では30%以内で測定値と一致した。Teschの式については、迷路入口の線量当量で規格化して比較したところ、全体として最大ファクター3程度で測定値と一致した。迷路入口におけるエネルギースペクトルの実験値と計算値の比較から、数MeV領域で計算値が過小評価となることを示した。熱中性子束分布に関しては、ターゲット室内同様に、迷路内でもファクター2程度計算値が過大評価となった。
小林 捷平*
JNC TJ9400 2000-009, 63 Pages, 2000/02
本研究の前半では、長半減期を有する核分裂生成物(FP)の代表的核種であるTc-99について、まず、中性子核データの現状、評価済核データについて調査を行った。次に、京都大学原子炉実験所の電子線型加速器と組み合わせて付設されている京都大学鉛スペクトロメータ(KULS)を用いて、99Tc(n,)100Tc反応断面積を熱中性子から1keVのエネルギー領域において測定した。中性子捕獲によって放出される即発ガンマ線は、アルゴンガス入り比例計数管を用いて測定した。入射中性子束/スペクトルの測定はBF3比例計数管によって行い、捕獲断面積の絶対値は10B(n,)反応によって求めた相対測定値をTc-99の熱中性子断面積に規格化した。Tc-99試料に対する中性子自己遮蔽効果は、MCNPコードを用いた計算によって補正した。Chouらが鉛スペクトロメータを用いて測定したTc-99の中性子捕獲断面積は、本測定値より全体に高いがエネルギー依存性はよく一致している。ENDF/B-VI、JENDL-3.2の評価済核データは、5.6、20eV共鳴領域及び数100keV以上で高くなっている。本研究の後半では、鉛スペクトロメータの原理とその構築、諸特性についてまとめた。中でも、KULSの特性として(1)中性子減速時間t(s)とエネルギーE(keV)の関係(ビスマス孔:E=190/tの2乗、鉛孔:E=156/tの2乗)及び(2)エネルギー分解能(ビスマス孔、鉛孔共に約40%)に関しては共鳴フィルターを用いて実験的に求め、(3)飛行時間分析法によってKULS体系中の中性子スペクトルを測定した。一方、MCNPコードを用いてKULSの諸特性を計算した結果、全体に実験値とよい一致を示した。
山本 敏久*; 北田 孝典*; 田川 明広; 丸山 学*; 竹田 敏一*
JNC TJ9400 2000-006, 272 Pages, 2000/02
多様な高速炉炉心の核特性に対する解析予測精度の向上を目的として、以下の3つの項目について検討を行った。第1部高速炉心の中性子スペクトルの誤差評価と計算精度向上策の検討高速実験炉「常陽」で用いられているスペクトルアンフォールディング法の精度を向上するため、初期推定スペクトル誤差を詳細に分析し、各々の誤差の大きさを定量的に評価するとともに、各誤差を積み上げることによって、より合理的な初期推定スペクトル誤差を評価することを試みた。検討の結果、初期推定スペクトル誤差に起因する誤差は相対的に小さく、断面積誤差に起因する誤差がほとんどであることがわかった。また、核分裂スペクトルの影響によって、数MeV以上の高速中性子束に無視できない量の誤差を生じることがわかった。第2部ガス冷却高速炉の解析手法に関する検討ガス冷却高速炉では、通常のNa冷却炉に比べて、冷却材チャンネルが体積割合に占める比率が大きく、顕著な中性子ストリーミング効果が現れることが予想される。一方、Na冷却炉用に提唱されている既存の手法では、冷却材チャンネルと平行な方向の拡散係数が無限大となり、そのまま適用することができない。本研究では、Kohlerが提唱した軸方向バックリングを考慮した方向依存拡散係数の概念を拡張し、ガス冷却炉でも正確に中性子ストリーミング効果が評価できる手法の検討を行った。第3部水冷却高速炉の解析手法に関する検討低減速の水冷却炉に対して、解析手法の違いによりどの程度計算結果に影響が現れるかについて検討を行った。軽水炉においては、燃料ペレット中の重核種の自己遮蔽効果が強い空間依存性を持つことが知られており、燃料ペレットを複数の領域に分割して評価する手法が用いられている。水冷却高速炉においても、冷却材として水を使用する以上、同様の問題が現れる可能性がある。検討の結果、燃料ペレット中の重核種の自己遮蔽効果の空間依存性は小さく、燃料領域を1領域として扱っても、臨界性、転換比ともに解析精度には問題が出ないことが確認された。
中根 佳弘
JAERI-Research 2000-003, p.110 - 0, 2000/02
中高エネルギー粒子場における放射線安全設計を行ううえで2次中性子による線量評価は重要である。本研究は数10MeV領域中性子入射における線量評価手法の確立を目的として、TIARAの40及び65MeV準単色中性子入射における平板ファントム内での中性子束及び線量を測定し、計算結果と比較した。また解析に必要な中性子線源スペクトル及び検出器応答関数を評価した。ファントム内中性子束及び吸収線量計算値は実験値をよく再現し、線量当量においても計算値は実験値をおおむね20%以内で再現したことから、このエネルギー領域における線量評価計算が20%程度の精度で可能であることを明らかにした。また固体飛跡検出器の応答では、開発した応答関数計算コード及び実験の結果から広範な中性子エネルギーに対する応答関数を決定した。
佐藤 聡; 飯田 浩正; Plenteda, R.*; Valenza, D.*; Santoro, R. T.*
Fusion Engineering and Design, 47(4), p.425 - 435, 2000/01
被引用回数:9 パーセンタイル:54.14(Nuclear Science & Technology)2次元及び3次元解析により、ITER/NBIポート周辺の遮蔽解析を行い、運転停止10秒後の生体線量率を評価した。運転中の2次元S解析、2次元放射化解析及び運転停止後の線に対する2次元S解析を行い、高速中性子束及び生体線量率分布を求めた。それらの値から、高速中性子束から生体線量率への変換係数を評価した。その結果、クライオスタット近傍では、1.5~410Sv/hour/(cmsec.)であることがわかった。トーラスの1/4を忠実にモデル化した。3次元モンテカルロ解析により、クライオスタット近傍の高速中性子束を求めた。分散低減技法の工夫により、統計誤差の小さい解が得られた。その結果、クライオスタット近傍の生体線量率は、20~100Sv/hourとなり、ITER/EDAの設計目標である100Sv/hourを満足することがわかった。
水野 正弘*; 宇都 成昭
JNC TN9400 98-007, 147 Pages, 1998/11
核燃料サイクル開発機構ではFBRの実用化に向けて必要となる各種炉内試験を実施するための、FBR安全性試験炉SERAPH(Safety Engineering Reactor for Accident PHenomenology)の設計研究を実施している。SERAPHでは定常及び種々の過渡試験が予定されており、試験時のドライバー炉心除熱のため重水が炉心冷却材に用いられる。種々の試験遂行には制御棒を使用するが、吸収体引抜き時のフォロワー材を重水とする制御棒案を候補の一つとして成立性の検討を行っている。この案では、重水の持つ高い減速比に起因して、重水フォロワー領域における中性子の減速・反射によって隣接燃料集合体で局所出力ピークが発生することが懸念されており、今後制御棒構造の具体化を進めるには詳細な局所出力分布特性評価手法の整備が不可欠である。このような背景に基づき、筆者らは制御棒周辺の局所出力分布を含む集合体内出力分布特性を適切に評価し得る核計算手法の整備を行った。解析ツールには制御棒近傍での中性子輸送挙動を統計的な影響を受けずに評価し得る2次元SN輸送計算コードTWOTRAN-IIを選定した。モデル化に際して、制御棒近傍での中性子平均自由行程を考慮し、制御棒とその周辺の15体分に相当するドライバー燃料集合体からなる2次元XY体系スーパーセルモデル、及びそれを13種類の単位セルモデルで構成する方法を考案した。スーパーセルモデルを効率良く構成するため領域マップ及びメッシュ境界を自動設定するプログラム、並びに計算で得られたメッシュ毎の中性子束をピン毎の出力密度に編集するプログラムを作成した。本成果は、今後制御棒構造やそれに関連した炉心構成の具体化に有効に活用されることが期待される。
野田 健治; Ehrlich, K.*; 實川 資朗; Moeslung, A.*; Zinkle, S.*
Journal of Nuclear Materials, 258(263), p.97 - 105, 1998/10
被引用回数:20 パーセンタイル:81.49(Materials Science, Multidisciplinary)国際核融合材料照射施設(IFMIF)の役割は、(1)核融合材料の開発、(2)核融合DESO炉の設計用材料データベースの取得等である。ここでは、IFMIFの概念設計活動(IFMIF-CDA)のためのユーザー要求を試験すべき材料、試験の形式、微小試験片技術、照射条件等について示される。試験すべき材料としては、フェライト鋼、V合金、SiC/SiC複合材料等の構造材料、LiO及びLiTiO等の増殖材料、炉内構造物絶縁及びRF窓等の絶縁セラミック材料等であり、試験の形式としては、照射後試験と運転条件下での材料特性を調べるための照射下その場実験がある。また、照射条件としては、必要な中性子束、照射試験体積、中性子束勾配、中性子スペクトル、材料損傷パラメータ等についてふれる。
鈴土 知明; Tuerkcan, E.*; H.Verhoef*
Nuclear Science and Engineering, 129(2), p.203 - 208, 1998/06
被引用回数:0 パーセンタイル:0.01(Nuclear Science & Technology)オランダボルセレ炉(PWR)の中性子束ゆらぎ解析を、ゆらぎの大きさに最も影響のある空間的にコヒーレントなゆらぎモードに注目して行った。解析の結果、0.03Hz付近に2次の方位角方向の振動モードの振動が観測された。また、0.1Hz付近にはコアワイドモードの振動が観測され、その強度は2次の方位角方向の振動モードに比べ非常に大きかった。冷却材温度等のほかの信号の解析もすることにより、2次の方位角方向の振動モードは熱水力ゆらぎから中性子束ゆらぎへの一方向的な影響によるもので、コアワイドモードは双方向、すなわち核熱カップリングによるものと示唆された。将来の大出力PWRでは冷却材温度ゆらぎの増大が懸念されるが、その研究結果は熱から核への一方向的な影響ではなく、双方向効果によるゆらぎモードを抑制することにより中性子束ゆらぎ量を抑えることができることを示唆している。
竹田 敏一*; 北田 孝典*; 山本 敏久*; 片木 洋介*
PNC TJ9605 98-001, 267 Pages, 1998/03
高速炉核特性の数値解析手法の改良として、マルチバンド法、摂動モンテカルロ法、輸送ノード法に関連する研究を行った。本報告書は以下の3部に分かれている。第1部 マルチバンド法による反応率計算法の改良マルチバンド法を用いて、ブランケット領域の反応率分布を詳細に評価する手法を検討した。フィッティング法によって作成した3バンドパラメータを用いて、U-238捕獲反応、U-235核分裂反応、Pu-239核分裂反応、U-238核分裂反応の反応率分布を解析した。対象核種としては、構造核種である鉄、ニッケル、クロム、およびナトリウムの4核種とした。マルチバンド法を用いることにより、いずれの反応率もブランケット深部で反応が増大する方向に補正され、補正量は最大で5%に達した。この結果は、従来の解析手法による実験値とのずれを改善する方向である。またこの補正量は、マルチバンド法におけるバンド間の散乱の取り扱いによって大幅に異なることがわかった。従来のフィッティング法の問題点を解決するべく、直接的なマルチバンドパラメータの作成法の検討も行った。第2部 摂動モンテカルロ法による反応度評価手法の改良摂動モンテカルロ法による摂動計算理論の検討及び、計算コードの作成を行った。昨年度までに使用していた相関サンプリング法だけでなく、導関数演算子サンプリング法でも計算できる、連続エネルギー摂動モンテカルロ計算コードを作成した。作成した計算コードを用いて「もんじゅ」炉心を対象とした計算を行い、参照解と比較検討した。「もんじゅ」にGEMまたは模擬燃料集合体を装荷した体系で、それらの集合体内のナトリウム密度を変化させた摂動、また制御棒全引き抜き体系で体系内のナトリウム密度を一様に変化させた摂動にともなう固有値の変化を調べた。ナトリウム密度の変化が小さい場合には、相関サンプリング法と導関数演算子サンプリング法のどちらの手法でも良好な結果を得ることができた。しかしながら、密度変化が大きい場合には、相関サンプリング法では妥当な結果を得ることができなかったが、導関数演算子サンプリング法では、そのような大きな密度摂動の場合でも良好な結果を得ることができることがわかった。第3部 3次元六角体系用輸送ノード法の改良集合体(ノード)内平均中性子束及びノード境界の中性子束から、集合体出力分布を評価する手法を、輸送理論に基づき導出し
Hunter
PNC TN9410 98-015, 81 Pages, 1998/02
本研究は、CEAから提供されたSuper-Phenixの起動試験炉心ベンチマークデータを動燃が解析した成果であり、動燃-CEA共同研究の一環として実施されたものである。動燃によるSuper-Phenixの解析結果を、CEAの解析結果及び実験測定値と比較したところ、CEAのC/E(解析/実験)値が系統的な径方向依存性を示すのに対して、動燃のC/E値はその3040%しかなく非常に小さいことが判明した。CEAが原因を検討した結果、両者のC/E値径方向依存性の違いの主たる要因は、使用した核データセット(JENDL-3.2CARNAVAL-IIII)にあると結論された。本検討の最終段階として、動燃はこの2種の核データセットの違い詳細に検討するために、感度解析を実施した。中性子束分布計算で用いた解析コードは2次元RZまたは3次元Hex-モデルのCITATIONとMOSESコードである。JENDL-3.2CARNAVAL-IIIIの違いに対する感度解析は、SAGEPコードを用いて行われた。ここでは、両者のエネルギー構造を統一するための縮約操作を施す必要があり、また、両者の核断面積の定義には幾つか食い違いがあることが分かった。感度解析の結果、JENDL-3.2とCARNAVAL-IIIIのC/E値径方向依存性の違いの原因は、少数の核種による寄与であることが判明した。両者の核データの比較結果は以下のとおりである。核分裂当たりの中性子発生数の違いは小さい(5%)。低エネルギーでの核分裂断面積差は大きい(30%、代表値10%)。下方散乱断面積は相対差としては大きい違いがあるが、絶対値の差は自群散乱と比較すれば無視できる。自群散乱の相対差は75%程度まであり、一般には20%以下である。捕獲断面積の違いは非常に大きく、30200%まで見られた。
山田 忠則; 相沢 雅夫; 白数 訓子
JAERI-Tech 97-072, 32 Pages, 1998/01
JRR-3Mのベリリウム反射体領域の垂直照射孔において、岩石型プルトニウム燃料試料の照射試験を原子炉運転4サイクル(平成7年1月23日~10月13日)に亘って行った。岩石型プルトニウム燃料試料は、トリア系化合物及びジルコニア系化合物をそれぞれディスク状に加圧成形し燃料ピンに収めてある。それぞれの燃料ピンは、キャプセル内に3段に分けて配置した。各段の目標試料温度を600C(上段)、800C(下段)、1000C(中段)となるように設計した。実際の照射試験では、原子炉出力20MWにおいて中段部の試料温度を1000Cに制御した。本報告書は、照射キャプセルの設計・製作及び照射試験の結果得られた試料の照射温度、中性子束に関する測定データをまとめたものである。
羽様 平
PNC TN9410 97-088, 139 Pages, 1997/10
重水臨界実験装置(DCA)の未臨界度測定試験体系の臨界特性(重水臨界水位、水位反応度、安全棒(制御棒)反応度価値、中性子束分布)を測定し、Sn法及びモンテカルロ法による計算結果と比較評価した。未臨界度測定試験体系は体系内にボイド領域を含むため、新型転換炉(ATR)の炉心核特性評価に用いられてきた拡散コードを適用することができない。ここではボイド領域を含む体系に適用できる輸送計算コードによる計算を行い、計算精度等の適用性を評価した。輸送コードとしては、2次元SN法によるTWODANTコードと、従来原子炉の計算にはあまり適用されていない多群モンテカルロ計算コードKENOを使用した。ATR等の原子炉炉心と比較して複雑な形状の未臨界度測定試験体系の炉心核特性を評価した結果、重水臨界水位等の臨界特性については、多群モンテカルロコードKENOの評価精度が0.5%k/k以内であることが確認され、中性子束の炉心内空間分布については、2次元SNコードTWODANTによって正しく評価できることが明らかとなった。特にモンテカルロコードは、炉心構造の違いによる微少な反応度変化を正確に再現しており、複雑な原子炉炉心の評価に極めて有効であることが確認された。
Fletcher, J. K.
PNC TN9410 97-065, 25 Pages, 1997/07
位置r、単位方向ベクトルの中性子束を(r, )と定義すると、多群輸送方程式は次式で表される。ただし、t(r)、はそれぞれ、全断面積、方向1 からへの散乱断面積、生成断面積を表し、また、は臨界係数を、(r,)は外部中性子源を表す。そして、この方程式を次の球面調和関数展開を用いて解く。ここで、(cos )はオーダーのルジャンドル陪関数で、とはそれぞれ方向の仰角及び方位角を表す。NはPN近似の次数を表す。三角関数の多項式である球面調和関数の直交性と漸化式を用いることにより、展開係数とに関する1階の微分方程式が導かれる。が奇数の項を消去することにより、拡散方程式の場合に用いられるような通常の有限差分法により解くことの可能な、2次の微分方程式が導かれる。メッシュ誤差低減は、その記述式の高次の差分項を保持したまま、2次式を用いて数値的に近似することにより行われる。当手法の採用により、メッシュ誤差は大幅に減少され、他の手法、特にモンテカルロ法により得られたものに匹敵する結果を直接計算することが可能となった。
竹田 敏一*; 木本 達也*; 北田 孝典*; 片木 洋介*
PNC TJ9605 97-001, 100 Pages, 1997/03
本報告書は次の2部と付録から構成されている。第1部 摂動モンテカルロ法による反応度評価手法の改良第2部 3次元六角体系用輸送ノード法の改良付録 高速炉におけるドップラー反応度解析のためのU238サンプルの実効断面積第1部 摂動モンテカルロ法による反応度評価手法の改良摂動モンテカルロ法の理論式の検討を行い、その後摂動モンテカルロ法の計算コードへの導入を行った。同じヒストリー数の計算を行ったところ、摂動モンテカルロ計算コードでの計算時間は、通常のモンテカルロ計算の12割程度の増加であった。作成した摂動モンテカルロ計算コードを用いて行った試計算結果は概ね妥当であり、また偏差も十分に小さいことから、摂動モンテカルロ法の有効性が示された。 しかしながら、得られた摂動前後の固有値の差が評価手法により、正や負になる場合があること、また、摂動による中性子源分布の変化を考慮しない従来手法と、摂動による中性子源分布の変化を考慮する新手法の間で、計算結果に有為な差が見られないことから、さらに摂動モンテカルロ計算コードに対して検討を加える必要がある。第2部 3次元六角体系用輸送ノード法の改良ノード法を用いた六角-Z体系用輸送計算コード「NSHEX」は、高速炉の炉心計算において非常に精度のよい評価を得ることがこれまでの研究で確かめられてきた。しかし非均質性の高い炉心においてややその精度が劣ることがわかっている。その原因として、ノード内空間分布を求める際用いる横方向もれの評価法が挙げられる。径方向スイープ時における、集合体からの径方向もれ分布を得るためには各ノード頂点中性子束を評価する必要がある。従来法では、その頂点の周囲のノード境界平均中性子束を用いている。新手法においては、その頂点近傍の中性子束分布を、いくつかの寄与が大きいと考えられるノード、およびノード境界の中性子束をパラメータとしてx,u 2次式で評価し、その分布式より頂点中性子束を算出している。以上の手法を用い、NEACRP 3D NEUTRON TRANSPORT BENCHMARKSの小型高速炉モデル、および実機「もんじゅ」体系を用いて検討計算を行った。その結果実効増倍率においては、多群モンテカルロ法によるGMVPに対して、どの手法もほぼ0.1%以内に一致する。各領域の中性子束も、数%以内に一致したものの、制御棒の挿入された体
吉田 昌宏; 長沖 吉弘
PNC TN9410 97-022, 34 Pages, 1997/02
「常陽」では、これまで、出力運転中の過剰反応度変化から燃焼係数(積算出力あたりの燃焼欠損反応度)を求め、解析結果との比較によりその予測精度を評価してきた。しかし、燃焼欠損反応度と核データあるいは中性子束計算精度との関係を詳細に評価するためには、上記方法で得られる炉心全体の特性ではなく、燃料集合体個々の燃焼反応度価値を測定・解析し、その予測精度と集合体の燃焼度や装荷位置の関係を明らかにする必要がある。このため、MK-II炉心を用いて燃焼反応度価値詳細測定試験を行い、個々の集合体の燃焼に伴う反応度の低下を測定した。試験では、燃焼度の異なる燃料集合体の置換反応度を測定することにより、集合体1体毎の燃焼反応度価値を評価した。燃料置換は燃焼初期(集合体平均燃焼度約1GWd/t)と寿命中期(約37GWd/t)及び燃焼初期と寿命末期(約62GWd/t)の2パターントし、測定位置は炉中心(第0列)、燃料領域中間(第2列)及び燃料/反射体境界(第4列)の3箇所とした。本試験により、集合体単位の燃焼反応度価値の燃焼度及び測定位置依存性を明らかにするとともに、運転監視コードシステム(MAGI)による計算結果との比較を行い、炉心管理上の燃焼係数予測精度を詳細に把握した。主な結果を以下に示す。(1)第0列の燃焼反応度測定結果は、燃焼初期から寿命中期までの燃焼で-0.19(%k/kk')、燃焼末期までの燃焼で-0.28(%k/kk')であった。(2) 燃焼反応度価値の径方向分布は、2つの燃焼パターンでほぼ一致し、第0列で規格化した時、第2列は0.67、第4列は0.28となった。(3)第0列と第2列の燃焼反応度価値の(MAGI)によるC/E値は良く一致し、燃焼初期から寿命中期までの燃焼で1.031.05、燃焼末期までの燃焼で0.940.95であった。また、反射体と隣接する第4列のC/E値は他2箇所と比較して57%高くなる傾向が見られた。 C/E値の燃焼度依存性に関しては、現在、測定に用いた燃料集合体のPIEを実施中であり、その結果に基づき詳細な評価を行う予定である。